Particle acceleration by circularly and elliptically polarised dispersive Alfven waves in 
a transversely inhomogeneous plasma in the inertial and kinetic regimes 
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Dispersive Alfven waves (DAWs) offer, an alternative to magnetic reconnection, opportunity to 
accelerate solar flare particles in order to alleviate the problem of delivering flare energy to denser 
parts of the solar atmosphere to match X-ray observations. Here we focus on the effect of DAW 
polarisation, left, right, circular and elliptical, in the different regimes inertial and kinetic, aiming to 
study these effects on the efficiency of particle acceleration. We use 2.5D particle-in-cell simulations 
to study how the particles are accelerated when DAW, triggered by a solar flare, propagates in 
the transversely inhomogeneous plasma that mimics solar coronal loop, (i) In the inertial regime, 
fraction of accelerated electrons (along the magnetic field), in the density gradient regions is ~ 20% 
by the time when DAW develops three wavelengths and is increasing to ~ 30% by the time when 
DAW develops thirteen wavelengths. In all considered cases ions are heated in the transverse to the 
magnetic field direction and fraction of heated particles is « 35%. (ii) The case of right circular, 
left and right elliptical polarisation DAWs, with the electric field in the non-ignorable transverse 
direction exceeding several times that of in the ignorable direction, produce more pronounced parallel 
electron beams (with larger maximal electron velocities) and transverse ion beams in the ignorable 
direction. In the inertial regime such polarisations yield the fraction of accelerated electrons 20%. 
In the kinetic regime this increases to 35%. (iii) The parallel electric field that is generated in 
the density inhomogeneity regions is independent of the electron-ion mass ratio and stays of the 
order O.Q3a;pe crrie/e which for solar flaring plasma parameters exceeds Dreicer electric field by eight 
orders of magnitude, (iv) Electron beam velocity has the phase velocity of the DAW. Thus electron 
acceleration is via Landau damping of DAWs. For the Alfven speeds of Va = 0.3c the considered 
mechanism can accelerate electrons to energies circa 20 keV. (v) The increase of mass ratio from 
mi/me — 16 to 73.44 increases fraction of accelerated electrons from 20% to 30 — 35% (depending 
on DAW polarisation). For the mass ratio mi/me — 1836 the fraction of accelerated electrons would 
be > 35%. (vi) DAWs generate significant density and temperature perturbations that are located 
in the density gradient regions. DAWs propagating in the transversely inhomogeneous plasma can 
effectively accelerate electrons along the magnetic field and heat ions across it. 

PACS numbers: 52.35.Hr; 52.35.Qz; 52.59.Bi; 52.59.Fii; 52.59.Dk; 41.75.Fr; 52.65. Rr 



I. INTRODUCTION 



Charged particles, accelerated to the speeds exceeding 
local thermal velocity, play a vital role in numerous lab- 
oratory, space, solar and astrophysical interdisciplinary 
processes. For example, in space plasmas, the Van Allen 
radiation belts form a torus of energetic charged plasma 
particles around Earth, which is held in place by Earth's 
magnetic field. These accelerated, charged particles come 
from cosmic rays. The radiation belts pose serious haz- 
ard to the space technology and future manned space 
travel. The outer electron radiation belt is mostly pro- 
duced by the inward radial diffusion [1[ and local ac- 
celeration y due to transfer of energy from whistler 
mode plasma waves to radiation belt electrons. Dis- 
persive Alfven waves (DAWs) are known to accelerate 
electrons in the Earth auroral zone [3]. In solar coronal 
plasmas, solar flare models involve magnetic reconnec- 
tion which converts the magnetic energy into heat and 
kinetic energy of accelerated particles. About 50-80% of 
the energy released during solar flares is converted into 
the energy of accelerated particles Q. Thus, knowledge 
of the details of particle acceleration via magnetic recon- 



nection (which is deemed a major mechanism operating 
solar and stellar flares) is of great importance. This is 
of direct relevance to the solar activity that affects hu- 
mankind via hazardous phenomena such as solar flares 
and energetic particles. In astrophysical plasmas, accel- 
erated particles play a major role in explaining the exis- 
tence of soft excess emission originating from clusters of 
galaxies that is of fundamental cosmological importance 
js!]. Also, multi- wavelength observations of Active Galac- 
tic Nuclei have now provided evidence that they are site 
of production of high temperature plasmas (kT > 100 
MeV) with strong tails of accelerated particles. The 
generation of hot, relativistic plasmas and accelerated 
particles is a challenging problem both in the context 
of Active Galactic Nuclei and the origin of high-energy 
cosmic rays [Q]. In laboratory plasmas, in the case of 
ions, plasma regimes in which enhanced fast-ion trans- 
port occurs because of strongly driven energetic particle 
models (EPMs) should be avoided in order to reach re- 
actor relevant performance and prevent first wall dam- 
age [7[ . The investigation of the transition from regimes 
characterised by weakly driven modes (such as Alfven 
Eigenmodes, e.g. toroidal Alfven eigenmodes, TAEs) to 



regimes dominated by strongly driven modes (like EPMs) 
and the determination of the related threshold are then 
of crucial importance for burning plasma operations [3] . 
In the case of electrons, e.g. in the laser driven inertial 
confinement fusion, accelerated electrons carry away a lot 
of energy pumped by the laser. Also, due to the laser's 
strong coupling with hot electrons, premature heating of 
the dense plasma (ions) is problematic and fusion yields 
are low. Significant improvements in the modelling of 
these phenomena have been achieved [a| . Edge Localised 
Modes (ELMs) are MHD instabilities destabilised by the 
pressure gradient in the H-mode Tokamak edge pedestal. 
ELMs produce losses of up to 10% plasma energy in sev- 
eral 100 micro-seconds. ELMs are major concern for 
the operation of ITER. Thus, ELM control is essential. 
ELMs are the focus of increasing attention by the edge 
physics community because of the potential impact that 
the large divertor heat pulses due to ELMs would have on 
the divertor design of future high power tokamaks such 
as ITER. Ref.[10| reviews what is known about ELMs, 
with an emphasis on their effect on the scrape-off layer 
and divertor plasmas. ELM effects have been measured 
in the ASDEX-U, C-Mod, COMPASS-D, DIILD, JET, 
JFT-2MJT-60U and TCV tokamaks, and were reported 
by Ref.[l3]- ELMs are known to produce super-thermal 
particle fluxes that are damaging for Tokamak wall ma- 
terials. 

At the beginning of the previous paragraph, interdis- 
ciplinary nature of the role of accelerated particles was 
emphasised. In all of the above mentioned fields of re- 
search, the mechanisms that lead to accelerated particles 
are in fact common, and can be split into three broad 
classes: (1) Electric field acceleration; (2) Shock acceler- 
ation and (3) Stochastic acceleration. 

In the electric field acceleration mechanism, Ref.|ll| 
considered dynamics of electrons under the action of two 
effects: the parallel (to background magnetic field) elec- 
tric field and friction between electrons and ions. The 
equation describing electron motion along the magnetic 
field is rrieVd = e_E — v^inieVd, where Vd is the electron 
drift velocity, Uei is the electron-ion collision frequency 
and dot denotes time derivative. In the limiting case of 
Vd << vthermai and the steady state regime, the equation 
of motion allows one to obtain the expression for Spitzer 
resistivity. In the limit of Vd > Vthermai, the steady 
state solution is not possible. This is because when the 
right hand side of the equation is positive i.e. when 
eE > VeiiTT-eVd^ One has electron acceleration. The latter 
leads to the so-called run-away regime. The acceleration 
due to the parallel electric field leads to an increase in the 
electron drift speed u^, in turn, this leads to a decrease in 
i^ei because i^ei oc vj ■ Therefore, when the electric field 
exceeds the critical value Ed = Uee^ In A/ {SiraQkBT) (so- 
called Dreicer electric field), the faster electron motion 
leads to the decrease in electron-ion friction, which in 
turn results in even faster motion. As a result the run- 
away regime is reached. For electric fields less than the 
Dreicer field, the rate of acceleration is less than the rate 



of collision losses and only a small fraction of the parti- 
cles can be accelerated into a non-thermal tail of energy 
E < LeEjj 0, where L is the acceleration lengthscale. 
Super-Dreicer fields, which seem to be present in many 
simulations of reconnection [12|, [ij] , accelerate particles 
at a rate that is faster than the collision or thermalisa- 
tion time. This can lead to a runaway and an unstable 
electron distribution which, as shown theoreticall y, by 
laboratory experiments and by the simulations [13, Il3| , 
will give rise to plasma waves or turbulence [ij, |l5| . 

In shock acceleration mechanism, which is also called 
first-order Fermi acceleration, the fractional momentum 
gain, Ap/p, per shock crossing is linear in the velocity 
ratio, Ush/v , where Ush and v are shock and particle 
speeds respectively. In stochastic acceleration mecha- 
nism, which is also called second-order Fermi accelera- 
tion, particles of velocity v moving along magnetic field 
lines with a pitch angle cos/x undergo random scatter- 
ing by moving agents with a velocity u. Because the 
head collisions (that are energy gaining) are more prob- 
able than trailing collisions (that are energy losing), on 
average, the particles gain momentum at a rate propor- 
tional to (m/w)^D^^, where D^^ is the pitch angle dif- 
fusion rate. Note that here dependence on the velocity 
ratio is quadratic (hence the name second order) [5|. 

In this work we focus on the electric field acceleration 
mechanism. In particular when the electric field is pro- 
vided by the DAWs. In magnetohydrodynamic (MHD) 
regime Alfven waves (AWs) cannot provide parallel elec- 
tric fields, therefore cannot accelerate particles via the 
electric field acceleration mechanism. After it was re- 
alised (la ) that in the kinetic regime (in which plasma is 
treated as a collection of charged particles, rather than 
a fiuid as in MHD) AWs with parallel electric fields are 
possible, a significant effort has been undertaken to study 
these waves. In laboratory plasmas, an externally ap- 
plied oscillating magnetic field (at a frequency near 1 
MHz for typical Tokamak parameters) was found to res- 
onantly mode convert to a kinetic AW, whose perpen- 
dicular wavelength is comparable to the ion gyroradius 
|17| . The kinetic AW, while it propagates into the higher- 
density side of the plasma after mode conversion, dis- 
sipates due to both linear and nonlinear processes and 
heats the plasma. If a magnetic field of 50 G effective 
amplitude is applied, approximately 10 MJ m~'^ of en- 
ergy can be deposited in 1 sec into the plasma. In space 
plasmas, surface waves, excited either by a MHD plasma 
instability or an externally applied impulse, were shown 
to convert in a resonant mode to a kinetic AW having a 
wavelength comparable to the ion gyroradius in the di- 
rection perpendicular to the magnetic field. The kinetic 
AW was found to have a component of its electric field in 
the direction of the ambient magnetic field and can ac- 
celerate plasma particles along the field line. A possible 
relation between this type of acceleration and the forma- 
tion of auroral arcs was discussed [I^. In solar coronal 
plasmas, it was found that an Alfven surface wave can be 
excited by 300 s period chromospheric oscillations which 



then resonantly excite kinetic AWs. Because of their ki- 
netic structure, the latter waves were found to dissipate 
their energy efficiently [l9l |. 

The parallel electric field is also produced when low fre- 
quency (w < UJc^^ where uJci = eB /rrii is the ion cyclotron 
frequency) dispersive AW has a wavelength perpendicu- 
lar to the background magnetic field comparable to any of 
the kinetic spatial scales such as: ion gyroradius at elec- 
tron temperature, Ps = y/k^iT^Jroi j tUci, ion thermal gy- 
roradius. Pi = ^JkBTi/mi/iOci^ |18j or to electron inertial 
length Ae = cjio^f^ [20|]- The central issue in understand- 
ing what physical effect (s) can provide the parallel elec- 
tric field can be understood based on the two-fluid theory 
of dispersive AWs 3] . The electron equation of motion, 

dve/dt+{ve-y)ve = ~{e/me)iE+VeXB)-V-Pe/{mene), 
indicates that in the linear regime, E'n can be supported 
either by the electron inertia term, dve/dt, or the elec- 
tron pressure tensor term, V • Pe/{mene). Dispersive 
Alfven waves are subdivided into Inertial AWs or Kinetic 
AWs depending on the relation between the plasma /3 and 
electron/ion mass ratio me/rrii. When /3 <^ me/rrii (i.e. 
when Alfven speed is much greater than electron and 
ion thermal speeds, Va ^ vth,i,vth,e) dominant mech- 
anism for sustaining E\\ is the parallel electron inertia 
(that is why such waves are called inertial AWs). When 
/3 ^ me/mi, (i.e. when Va <C vth,i,vth,e) then clearly 
thermal effects become important. Thus, the dominant 
mechanism for sustaining £'|| is the parallel electron pres- 
sure gradient (that is why such waves are called kinetic 
AWs - kinetic motion of electrons is the source of pres- 
sure). In the MHD regime, it has been known [21j that 
the phase-mixing mechanism, in which an AW moves 
along the density inhomogeneity that is transverse to 
the background magnetic field, can provide a fast dis- 
sipation mechanism. Fast in a sense that AW ampli- 
tude decrease, as a function of travel distance, x, along 
wave propagation direction, is a much faster function 
of X, B±{x) — exp[(— 77V'^^fc^/6V'^)a;'^], compared to the 
usual resistive damping law Bj_{x) — exp[(— Tyfc^/VA)^], 
where 77 is the plasma resistivity and V^ denotes spa- 
tial derivative of the Alfven speed in the transverse to 
the background magnetic field direction 21]. The fast 
dissipation occurs because the phase mixing effect leads 
to the build-up of sharp gradients in the Alfven wave 
front, in the transverse to its propagation direction. It 
has been recently discovered [22] in the solar coronal con- 
text that when the phase-mixing mechanism is consid- 
ered in the kinetic regime, electrons can be effectively 
accelerated by the phase-mixed dispersive AW. It was 
found 22] that when dispersive AW moves in the trans- 
versely inhomogeneous plasma, parallel electric field is 
generated in the regions of plasma inhomogeneity. This 
electric field can effectively accelerate electrons for a wide 
range of plasma parameters [23] . It turns out that simi- 
lar electron acceleration mechanism has been also found 
in the Earth magnetosphere context [2J, [23|. See also 
Ref . [26l for the cross-comparison. The difference is that 
Ref. [22,[2J| have studied the case of plasma over-density. 



transverse to the background magnetic field, whereas 
Ref. j24| - |26t considered the case of plasma under-density 
(a cavity), as dictated by the different applications con- 
sidered (solar coronal loops and Earth auroral plasma 
cavities). Ref.[27| studied the kinetic-scale Alfven wave 
phase mixing at the boundaries of the density cavity, the 
strongly varying Alfven speed profile above the Earth au- 
roral ionosphere (so-called ionospheric Alfven resonator) , 
which leads to the development of the parallel electric 
fields that accelerate electrons in the aurora. They per- 
formed simulations to study the evolution of these fields 
including both parallel and perpendicular inhomogeneity. 
Calculations of Ref. [28|] , based on the linear kinetic the- 
ory of Alfven waves, also indicate that Landau damping 
of these waves can efficiently convert the Poynting fiux 
into field-aligned acceleration of electrons. More recently, 
an analytical theory has been developed I29ll to explain 
the numerical simulations presented in Ref. |23. [23| . Also, 
a possibility of acceleration of electrons by inertial AWs 
has been explored in a one dimensional, homogeneous 
plasma case (without phase mixing) [30|. The latter ex- 
tended and adopted for solar coronal context earlier work 
of Ref. [31| who demonstrated the Fermi-like acceleration 
by inertial Alfven waves of a small proportion of elec- 
trons, up to speeds of twice the Alfven speed and argued 
that this is a plausible scenario for auroral electron accel- 
eration. The effect of the Hall term on the phase mixing 
of AWs has been studied in Ref. [S^l . 

This work aims to fill the gap in understanding of 
particle acceleration by dispersive Alfven waves in the 
transversely inhomogeneous plasma via full kinetic sim- 
ulation particularly focusing on the effect of polarisation 
of the waves and different regimes (inertial of kinetic). 
In particular, we study particle acceleration by the low 
frequency (a; = O.Swci) DAWs, similar to considered in 
Ref. [23, [23 , but here we focus on the effect of the wave 
polarisation, left- (L-) and right- (R-) circular and ellip- 
tical, in the different regimes inertial (/? < mp^/mi) and 
kinetic (/3 > me/mi). The primary goal is to study the 
latter effects on the efficiency of particle acceleration and 
the parallel electric field generation. 



II. THE MODEL 

We use EPOCH (Extendible Open PIC Collaboration) 
a multi-dimensional, fully electromagnetic, relativistic 
particle-in-cell code which was developed and is used 
by Engineering and Physical Sciences Research Council 
(EPSRC)-funded Collaborative Computational Plasma 
Physics (CCPP) consortium of 30 UK researchers. We 
use 2.5D version of the EPOCH code which means that 
we have two spatial component along x- and y-axis and 
there are three Vx,Vy,Vz particle velocity components 
present (for electrons and ions). The relativistic equa- 
tions of motion are solved for each individual plasma 
particle. The code also solves Maxwell's equations, with 
self-consistent currents, using the full component set of 



^ in the least dense parts of the domain 
= Vmax), which are located at the lower- 



Here uj. 



EM fields E^,Ey,E^ and B^,By,B:,. EPOCH uses un- 
normalised SI units, however, in order for our results 
to be general, we use the normalisation for the graph- 
ical presentation of the results as follows. Distance and 
time are normalised to c/ujpe and cjpg^ , while electric and 
magnetic fields to Upecme/e and WpeWe/e respectively. 
Note that when visualising the normalised results we use 
no = lO^'^ m" 
(2/ = and y 

most and uppermost edges of the simulation domain (i.e. 
fix w„p — 5.64 X 10^ Hz radian on the domain edges). 

= •\/«-ee^/(eom.e) is the electron plasma fre- 
quency, Ua is the number density of species a and all 
other symbols have their usual meaning. The spatial di- 
mension of the simulation box are varied as a; = 5000 
and y — 200 grid points for the mass ratio mi/nif. — 16, 
and X — 10712 and y = 200 grid points for the mass 
ratio mi/me = 73.44 which represents l/25th of the 
actual irii/me of 1836. The choice of the mass ratio 
mi/nie = 1836/25 — 73.44 is motivated by following rea- 
sons: (i) since the DAWs considered in this study are 
essentially Alfvenic (for the fixed throughout this paper 
value oi UJ = OSaJd there is very little difference between 
Alfven branch with uj = VAk and L- and R- circularly 
polarised DAW branches that propagate along the mag- 
netic field, see e.g. Fig. 3.4 from Ref.fSJ]), as the Alfven 

1/2 

seed Va = Bo/y/JI^yiftpn scales as m^ the actual Alfven 

wave would be 5 times slower than considered here; (ii) 
rrii/me = 1836/25 = 73.44 value lands us in the kinetic 
regime because plasma beta in this study is fixed at /3 = 
2{vth,^/cf{wp,/uJ,,f = noiO,0)kBT/{B^/{2t^o)) = 0.02. 
Thus 13 = 0.02 > me/mi = 1/73.44 = 0.0136; (iii) This 
is the maximal value that can be considered with the 
available computational resources - using 32 Dual Quad- 
core Xeon = 32 X 8 = 256 processor cores. A typical 
with mi/me — 73.44 run takes 48 hours on 256 proces- 
sors. The grid unit size is A^i. Here A^i ~ vtt.e/'^pe is 
the Debye length {vth,e = \/kBT jra^ is electron ther- 
mal speed). We use 100 electrons and 100 ions per cell 
which means that for nii/me = 16 there are 2 x 10* and 
for mi/nie = 73.44 there are 4.2848 x 10* particles in 
the simulation box. Plasma number density n{x, y) (and 
hence ujpe{x,y)) can be regarded as arbitrary, because 
we wish our results to stay general. We impose con- 
stant background magnetic field Bqx — 320.753 Gauss 
along X-axis. This corresponds to Bq^ — 1.0(wpem-e/e), 
so that with the normalisation used for the visualisa- 
tion purposes, normalised background magnetic field is 
unity. This fixes ujce/ujpe — 1.00. Electron and ion 
temperature at the simulation box edge is also fixed at 
r(0,0) = Te(0,0) = T,(0,0) = 6 X lO^K. This in con- 
junction with no (0,0) ~ 10^^ m^'^ makes plasma pa- 
rameters similar to that of a dense flaring loops in the 
solar corona. For example radiative hydrodynamic mod- 
elling of the Bastille-Day solar flare by Ref. (3^ (see their 
Fig. lb and Ic) indicate loop apex temperatures and den- 
sities as 5.3 X lO^K and 2.5 x 10^"^ m~^ respectively. This 
is consonant with our overall goal to provide a viable. 



TABLE I: Numerical simulation parameters 



Regime 


Inertial 


Kinetic 


mi/mi, 


16 


73.44 


Uce/^pe 


1.000 


1.000 


P 


0.020 


0.020 


C/U)pe [m] 


0.053 


0.053 


Ad = J'i,e[m] 


0.005 


0.005 


Vth,e/C 


0.101 


0.101 


Vth,i/C 


0.025 


0.012 


Va/c = bJci/uJpi 


0.25 


0.117 


Va.vh/c 


0.243 


0.116 


Vl/c 


0.201 


0.097 


Vr/c 


0.264 


0.131 


tend = 75a;;/ [x IQ-^s] 


2.127 


9.763 


Uy 


200 


200 


rix 


5000 


10712 



fully-kinetic model for particle acceleration by means of 
DAWs during solar flares to alleviate so called "number- 
problem" . The latter refers to the too high total number 
(10'^'* — 10^^ per second) of accelerated electrons required 
to produce the observed hard X-ray emission compared 
to that available in the corona, if the particle acceleration 
takes place at the loop apex. This would mean that if the 
solar flare particle acceleration volume is in the range of 
1 — 10 Mm'^ with the number density of n = 10^^ m~^, to 
match the observational lO'^* — 10'^^ accelerated electrons 
per second, full 100% of electrons need to be accelerated. 
There is no known particle acceleration mechanism that 
operates with the 100% efflciency. Thus it would seem 
attractive to have flare energy release triggering DAWs 
which rush down the footpoints and accelerate electrons 
along the way to chromosphere [35|. Such scenario is 
rather attractive in the light of the earlier work Ref. [22 1 
which showed that such particle acceleration is indeed 
possible. The summary of numerical simulation parame- 
ters is given in Table ID 

As in Ref.[22|, in order to model solar coronal loops, 
we consider a transverse to the background magnetic field 
variation of number density as following 



ni{y) = Ueiy) = H-3exp 



y - lOOAz; 

50An 



^ f{y)- 

(1) 

Eq.(l) implies that in the central region (across the y- 
direction), the density is smoothly enhanced by a factor 
of 4, and there are the strongest density gradients having 
a width of about 30Ad around the points y — 51.5A£i 
and y = 148. 5Ad. The background temperature of ions 
and electrons, are varied accordingly 



T,{y)/To = TM/To = fiyy 



(2) 



such that the thermal pressure remains constant. Be- 
cause the background magnetic field along the x- 
coordinate is constant, the total pressure is constant too. 
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FIG. 1: Contour (intensity) plots of the following physical quantities at tend = 75a;~^ (the time when DAW develops three 
wavelengths): (a) Bx{x,y,t = tend) - -Bo, (b) By{x,y,t = tend), (c) B^{x,y,t = tend), (d) E^{x,y,t = tend), (e) Ey{x,y,t = 
tend), (f) E^{x,y,t = tend), (g) Jx{x,y,t = tend), (h) Jy{x,y,t = tend) and (i) J^{x,y,t = tend)- See text for the normalisation 
used. This figure pertains to the numerical run L16 (see Table lU for the run parameters). 



Thus the pressure balance in the initial conditions is fully 
enforced. Further we drive domain left edge, x — IXd, 
as follows 



Eyil,y,t + At)^Eyil,y,t)- 
Ay sm{ujdt) (1 - exp [-(i/io)^]) , 

E,il,y,t + At)=E,il,y,t)- 
A, cos{ujdt) (1 - exp [-{t/to)^]) , 

for the case of left-polarised DAW and 

Ey{l,y,t + At)=Ey{l,y,t) + 
Aysm{ujdt) (1 - exp [-(^Ao)^]) , 



E,{l,y,t + At) = E,il,y,t)~ 
A, cos{ujdt) (1 - exp [-{t/ta)^]) , 



(3) 



(4) 



(5) 



(6) 



for the case of right-polarised DAW. Here uJd is the driv- 
ing frequency which was fixed at w^ = O.Swci (Although 
numerical value of ujd actually varies when ion mass is 
varied). This ensures that no significant ion-cyclotron 
damping is present and also that the generated DAW is 
well on the Alfven wave branch with dispersion proper- 
ties similar to the Alfven wave (see e.g. Fig. 3.4 from 



Ref . j33| ) . to is the onset time of the driver, which was 
fixed at 3.0uj~^ i.e. 48.000ajZ,^ for the case oimi/nie — 16 
and 220.320ujp^^ for the case of mi/rrie = 73.44. This 
means that the driver onset time is about 3^ = 9w^^. 
Imposing such a transverse electric field E± = Ey^z driv- 
ing on the system results in the generation of L- or R- 
circularly or elliptically polarised DAW. Definition for 
the L- or R- polarisation is a follows: Consider the mag- 
netic field directed along positive x-axis direction and 
electric field vector E rotates in j/z-plane. In the case of 
L-polarisation, if left hand thumb points along the mag- 
netic field (and hence along the propagation direction of 
DAW) curled fingers show the direction of rotation of E, 
i.e. clockwise. In the case of R-polarisation then right 
hand thumb points along the magnetic field and curled 
fingers show the direction of rotation of E, i.e. counter- 
clockwise. Note that ions rotate in the same direction as 
E of the L-polarised DAW (that is why ion-cyclotron res- 
onance can occur), at the same time electrons rotate in 
the opposite sense. For the DAW the dispersion relation 
is as following (see Eqs.(3.78) and (3.87) from Ref. 33]): 



1/2 



fc=^ 



pe 



pi 



(Uce ±i^){uJci T^) 



(7) 



where upper signs are for the L-polarisation and lower 



signs for the R-polarisation. Eq.(7) shows that for L- 
polarisation the allowed frequency range of < w < Wci 
which indicates existence in the ion-cyclotron resonance 
at w = ujci- Ions can resonate because they rotate in the 
same direction as the electric field vector of the DAW. 
For the R-polarisation the allowed range is < cj < LOce 
with the electron-cyclotron resonance aX uj — ujce- In the 
plasma literature L-polarised wave is usually referred to 
as ion cyclotron wave because it cuts off at uj — ujd, while 
R-polarised wave is referred to as the whistler branch 
as it cuts off at w = uJce- We follow the terminology 
adopted in [3| and refer to all waves with perpendicular 
wavelengths comparable to any of the above mentioned 
kinetic scales as DAWs, irrespective whether they are 
on L- or R-polarisation dispersion branch and inertial or 
kinetic regime. For uj < O.Sujci from the dispersive prop- 
erties point of view there is little difference between L- 
and R-polarised waves and the shear Alfven wave which 
has the phase speed of VA,ph/c = iyAJc)! ^X^ViyAJc^ 
which for the parameters considered is slightly different 
from the Alfven speed VaJc = ujd/ujpi- Eq.(7) has been 
used to calculate L- and R- polarisation wave phase seeds 
in Table n 



n, (a) 



ni(b) 
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e,5e+15 2.1e+16 3,6e+16 5.0e+16 e,7e+15 2.1e+16 3,5e+16 5.0e+16 
T, (c) T; (d) 



1D0 200 3D0 400 500 100 200 300 400 500 
» [c/"„,] » [c/"„,] 
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FIG. 2: Contour (intensity) plots of the following physical 
quantities at tend = 75ct;~^ (the time when DAW develops 
three wavelengths): (a) ne{x,y,t = tend), (b) ni{x,y,t = 
tend), (c) Te{x,y,t = tend), and (d) Tt(x,y,t = tend)- Here 
physical quantities are in SI units. This figure pertains to the 
numerical run L16. 



The initial amplitudes of the E± are chosen as follow- 
ing: (i) in the case of circular polarisation Ay = Az — 
0.05(cB^o); Thus, the relative DAW amplitude is 5% 
of the background, thus the simulation is weakly non- 
linear, (ii) in the case of elliptical polarisation when 
electric field in the ignorable z direction is six times 
larger than in y direction, Ay = O.Ob^cBxo) x ^^2/37 
and Az = 0.05{cB^o) x ^^2/37 x 6; (in) in the case of 
elliptical polarisation when electric field in the ignorable 
z direction is six times smaller than in y direction, Ay — 
0.05{cB^q) X ^2/37 X 6 and A^ = 0.05{cB^o) x ^/2J37. 



TABLE II: Numerical simulation run identification and phys- 
ical parameters. I stands for inertial and K for kinetic. Reg. 
stands for regime. 



Run ID 


Polaris. 


mi/nie 


Reg. 


Ay/A-. 


iendKl^] 


Figs. 


L16 


L — circ. 


16 




1 


75 


1,2,3 


R16 


R — circ. 


16 




1 


75 


4 


EL16 


L - ellip. 


16 




6 


75 


5 


ER16 


R - ellip. 


16 




6 


75 


6 


ELI61 


L — ellip. 


16 




1/6 


75 


7 


ERI61 


R - ellip. 


16 




1/6 


75 


8 


L73 


L — circ. 


73.44 


K 


1 


75 


9 


R73 


R — circ. 


73.44 


K 


1 


75 


10,11,12 


EL73 


L — ellip. 


73.44 


K 


6 


75 


13 


ER73 


R - ellip. 


73.44 


K 


6 


75 


14 


EL73i 


L - ellip. 


73.44 


K 


1/6 


75 


15 


ER73i 


R - ellip. 


73.44 


K 


1/6 


75 


16 


L16Long 


L — circ. 


16 


I 


1 


300 


17,18 



Note that the amplitudes Ay and Az have been chosen 
such that in all considered cases of circular or ellipti- 
cal polarisation, the wave power pumped into the system 
is the same; i.e. in all presented runs Ay + A^ = const, 
with the constant being the same for each numerical run. 
This ensures consistency when studying the efhciency of 
particle acceleration and i?|| generation. 

The considered numerical runs with their identifying 
name used throughout this paper is shown in Table HIl 



III. RESULTS 

Below we present numerical simulation results for the 
runs given in Table HIl Run L16 is similar to one consid- 
ered in |22|. The latter used numerical code TRISTAN. 
Throughout this paper we use numerical code EPOCH. 
The run L16 provides important benchmarking and vali- 
dation for our results. We gather from Fig. 1(e) and 1(f) 
that transverse electric field driving, as described above, 
generates DAWs that for < a: < 250 propagate to the 
right, i.e. along the magnetic field. However, due to the 
driving at x = 1 also generates a wave that propagates to 
the left, because of the periodic boundary conditions used 
both for particles and fields throughout this paper, the 
left-propagating wave appears on the right hand edge of 
the domain and propagates in the domain 250 < x < 500. 
In turn, such driving also generates transverse magnetic 
fields which are shown in Figs. 1(b) and 1(c). One can 
clearly see from Figs. 1(c) and 1(e) that DAW compo- 
nents {Ey,Bz) experience phase mixing. That is initially 
flat wave front becomes progressively distorted, because 
in the middle of the domain phase speed of the wave is 
slower than at the edges. Figs. 1(b) and 1(c) show that 
DAW has developed three full wavelength. Period of the 
DAW is T = 27r/a;d = 27r/(0.3Lj„) « 21uj-^. Three pe- 
riods (and hence three wavelengths) will take 63aj~^ to 
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FIG. 3: Time evolution (at f = and t = te„d = 75a;^^) of electron and ion velocity distribution functions (on a log-linear 
plot): (a) fe{vx,t = 0) dashed (inner) curve and fe{vx,t = tend) solid (outer) curve, (b) fe{vy.,t — 0) dashed (inner) curve and 
fe{vy, t = tend) solid (outer) curve (note that the both curves overlap), (c) fe{vz,t = 0) dashed (inner) curve and fe{vz,t = tend) 
solid (outer) curve (note that the both curves overlap), (d) fi{vx,t — 0) dashed (inner) curve and fi{vx,t = tend) solid (outer) 
curve (note that the both curves overlap), (e) fi{vy,t — 0) dashed (inner) curve and fi{vy,t = tend) solid (outer) curve, (f) 
fi{vz,t = 0) dashed (inner) curve and fi{vz,t = tend) solid (outer) curve. Time evolution (at 20 time intervals between t = 
and t = tend) of the following: (g) ma.x\Ex{x,y,t)\, triangles connected with a solid curve and (h) AH^^^eit) index, diamonds 
connected with dashed curve, according to Eq.(8), AHj_^i{t) index, triangles connected with a solid curve, according to Eq.(9). 
This figure pertains to the numerical run L16. 



develop which added to the DAW driver on-set time of 
9u}~^ gives 72(jj~^ which roughly coincides with the end of 
the simulation time of 75aj^^ (shortfall of 75 — 72 = 3aj^^ 
could be roughly attributed for the latency time it takes 
for ions to respond to DAW driving). 

We gather from Figs. 1(a) and 1(d) that in the regions 
where density has transverse gradients longitudinal (par- 
allel) magnetic field (B^) and electric field (£'|| = E^) are 
generated. It can be clearly seen that these are generated 
only on the density gradients located within 4 < y < 7 
and 13 < y < 16. Outside these ranges as can be corrob- 
orated by plotting n{y) from Eq.(l) the plasma density 
is constant. This means the transverse gradient scale is 
7 — 4 = 16 — 13 = 3c/wpe ~ 30Xd for mi/rrie = 16 corre- 
sponds to 0.75c/ujpi, i.e. the transverse density gradient 
is about one ion-inertial lengths long. Further insight 
can be obtained by looking at x, y, z components of the 
plasma current J = e{niVi — UeVe). Our normalisation 
for the current is Jnorm = &{ne){vth.e) where (rig) is the 
number density averaged over the transverse, y-direction 
and {vth,e) = \/kB{Te)/me with (Tg) being electron tem- 
perature also averaged over y-direction. In effect Jnorm 
is a typical electron thermal velocity current. We gather 
from Fig. 1(g) that sizable longitudinal currents (« 0.8) 
are generated in the transverse density gradient regions. 
These currents provide explanation for the source of the 



i?l|. As also outlined by Ref.jSg], charge separation be- 
tween more mobile electrons and less mobile ions causes 
generation of time variable i?|| . The transverse currents 
Jy (Fig. 1(h)) and Jz (Fig. l(i)) are manifestations of 
the driving DAW which actually follow the Ey and E^ 
spatial pattern respectively. 

Fig. 2 displays behaviour of number density and tem- 
perature of electrons and ions. The noteworthy features 
are: (i) the transverse (kink-like) oscillations of plasma 
density that are mainly confined the density gradient re- 
gions (Figs. 2(a) and 2(b)); (ii) hot plumes of plasma 
w 8.7 X lO^K that emanate from the density (and tem- 
perature) gradient regions where particle acceleration by 
DAWs takes place (Figs. 2(c) and 2(d)). Fig. 2 also gives 
a good indication how the initial pressure balance looks 
like: number density goes up from 10^^ cm~'^ on the do- 
main top (y = Ly^max = 20.12) and bottom (y = 0) 
edges and increases to 4 x 10^® cm~'^ in the middle of the 
y-range. 

Fig. 3 shows the distribution function time evolution, 
quantifies the generated parallel electric field amplitude 
and associated particle acceleration for the run L16. We 
gather from Figs. 3(a)-3(f) that as in [24I the gener- 
ated parallel electric field in the density in homogeneity 
regions produces electron beams propagating along the 
magnetic field (Fig. 3(a)), while in the transverse direc- 
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FIG. 4: As in Fig. 3 but for the numerical run R16. 
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FIG. 5: As in Fig. 3 but for the numerical run EL16. 



tions y and z the electron distribution function remains 
unchanged (Fig. 3(b) and 3(c)). In the case of ions only 
temperature increase is seen in the transverse direction to 
the magnetic field. Note that Ref.[33| also found heating 
of ions with L-circularly polarised AWs using test-particle 
simulations in the context of solar chromospheric heating. 
It should be noted that Ref.[37] considers non-resonant 
wave-particle interactions, in which the polarization, as 
such, is not important. No ion beam have developed 
as we see only broadening of the distribution function 



(Fig. 3(e) and 3(f)). Along the magnetic field the ion 
distribution function remains unchanged (Fig. 3(d)). The 
evolution of the amplitude of the parallel electric field is 
studied in Fig. 3(g) where we plot max \Ex{x, y,t)\. Since 
parallel electric field is mostly confined to the density 
gradient regions Taeoi\Ex{x,y,t)\ as function of time es- 
sentially tracks the amplitude of the generated parallel 
electric field in the density gradients. We see that after 
initial growth it saturates at about 0.03. We quantify the 
particle acceleration and plasma heating by defining the 



following indexes: 



^^||,e(i) 



^I^IXt-th,.) f'^i'"^^*)^'"^ / i'^LjH,y X Lr,^,nax) 



J — QQ Je\^X7^)^^x I \-'-^y,rnax -^ ^x^max) 



(8) 



AH^,,{t) 



/k^k^,,) Mv^,t)dv^ / {2LiH^y X L^, 



IToo fii'"-Li^)d'V± / {Ly^max X Lx^max) 



(9) 
Here fe,i are electron or ion velocity distribution func- 
tions. They are normalised in such a way that when in- 
tegrated by all spatial and velocity components they give 
the total number of real electrons and ions (not pseudo- 
particles that usually mimic much larger number of real 
particles in the Particle- In-Cell simulation). LjH,y is the 
width of each density gradients according to Eq.(l), lo- 
cated within 4 < y < 7 and 13 < y < 16. Ly^rnax is the 
full width in the y-direction. The numerator of Eq.(8) 
gives the number of superthermal electrons in the paral- 
lel to the magnetic field direction, as a function of time, 
divided by the area (2LiH,y x Lx^max) where the electron 
acceleration takes place. As we know from Fig. 1 parallel 
electric field is generated only in the two density gradi- 
ent regions. The denominator gives total number of real 
electrons in the whole simulation box, at i = 0, divided 
by the total area (Ly^max x L^jnax)- Thus, the defini- 
tion used in Eq.(8) effectively provides the fraction (the 
percentage) of accelerated electrons in the density gradi- 
ent regions. The same logic applies to ions as refiected 
by Eq.(9) but in the transverse direction. The choice 
for such indexes has several reasons: (i) AH^^_f,{t) and 
AHj_^i(t) quantify fraction of electrons and ions at time 
t with velocities greater than average thermal velocity 
in density gradient regions, in the parallel and perpen- 
dicular direction respectively; (ii) As can be seen from 
Figs. 3-9 and 12-16, driving by DAWs mostly affect elec- 
tron motion parallel to the magnetic field, while for ions 
the effect is only in the perpendicular direction; (iii) They 
way how A7J|| e(t) and AH±^i{t) are defined alleviates 
the problem that the distribution function value varies 
as the total number of particles in the simulation is al- 
tered in different runs (because domain size needs to be 
adjusted so that DAWs do not collide producing interfer- 
ence). For example comparing Fig. 3(a) and 9(a) we see 
that the height of the peak at Vx.e = is different despite 
the fact that both graphs show the distribution function 
for the medium with the same temperature distribution 
and other parameters. This is due to the fact that distri- 
bution function is normalised such that when integrated 
by all spatial and velocity components, it gives the total 
number of particles (not pseudo-particles) in the simula- 
tion. Fig. 9 pertains to the run L73 in which number of 



grids {ux = 10712) is more than twice compared to the 
run L16 in which Hx — 5000. Since we always load 100 
electrons and 100 ions per cell, this results in doubling the 
peak height in Fig. 9(a) compared to Fig. 3(a), as there 
are twice as many particles in run L73 than in run L16. 
It is clear from the definition (Eq.(8)) that this difference 
in heights cancels out, thus providing the true precent- 
age of accelerated electrons in the density gradient re- 
gions. The only drawback our Aif|| e(t) and AH±^i{t) in- 
dexes have is that they cannot distinguish between parti- 
cle acceleration (in which case bumps on the distribution 
function appear) and heating (in which case distribution 
function simple broadens, proportionately retaining its 
shape at t = 0). Despite this AiJ|| g(i) and AH±,i{t) 
indexes still provide useful information about the frac- 
tion of particles at time t with velocities greater than 
average thermal velocity in the density gradient regions. 
We gather from Fig. 3(h) that AiJ||.e(t) and AH±^,{t) 
indexes attain values 0.2 and 0.35, respectively. This 
means that in the density gradient regions 20% of elec- 
trons are accelerated in the parallel to the magnetic field 
direction, and 35% of ions gained velocity above thermal 
(were heated up) in the transverse direction. It should 
be noted that these conclusions apply for the end simula- 
tion time oitend — '^^^ci'- See discussion of run LlGLong 
which discusses long term dynamics that corresponds to 
tend = SOOwj/. 

It is interesting to note that in Ref. [23| we used much 
more computationally demanding definition (see Eq.(5) 
from Ref. [2 3]). The demanding part is storing large 
data for the distribution function {J{vx^e,x,y)) in order 
to perform integration by dxdy. Our above definition 
(Eq.(8)) is much less demanding because it only stores 
f(vx,e) and spatial integration is replaced by the factor 
l/{2LiH.y/Ly^rnax)- By Comparing the rightmost data 
point for AH^\,,{t) in Fi g. 3 (h) with middle data point 
from Fig. 2(b) from Ref. [23| (because their physical pa- 
rameters are the same), we gather that they both yield 
value of 0.2, thus proving that the both definitions are 
quivalent. Another observation we make is that AiJ|| e(t) 
starts from about 0.15. Which is reasonable because it 
is known that for a Maxwellian distribution (which holds 
before the particle acceleration takes place) 16% of par- 
ticles have velocty above thermal. This, of course, ap- 
plies to a Maxwellian with a uniform temperature. In 
our case plasma temperature varies across y coordinate, 
thus the Maxwellian does not appear as a parabola in 
e.g. Fig. 3(a), as is it should have done because it is a 
log-linear plot. 

Figs. 4, 5 and 6 correspond to the runs R16, EL16 and 
ER16. We group them together because their behaviour 
is similar, but notably different from run L16 (Fig. 3). 
We see from Figs. 4(a), 5(a) and 6(a) that bumps on the 
electron parallel velocity distribution function are much 
more pronounced and they now extend above 0.4c. 

In the case of ions we gather from Figs. 4(f), 5(f) and 
6(f) that transverse beams start to form in the ignorable z 
spatial direction. As far as the generated parallel electric 
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FIG. 6: As in Fig. 3 but for the numerical run ER16. 
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FIG. 7: As in Fig. 3 but for the numerical run ELI61. 



field, in the density gradient regions, it saturates at about 
the same level of 0.03. yl7J|| e(i) and AH±^i{t) indexes 
also exhibit similar behaviour. This is due to the fact 
that more efficient electron and ion acceleration produces 
relatively small number of more energetic particles, which 
do not affect in the indexes greatly. 

Intermediate conclusions that can be drawn from com- 
paring runs L16, R16, EL 16 and ER16 are that in the in- 
ertial regime, R-circularly polarised, L-and R-elliptically 
polarised DAW with electric field in the non-ignorable 



(y) transverse direction exceeding 6 times that of in the 
ignorable (z) direction produce more pronounced parallel 
electron beams (in v^) and transverse ion beams in the 
ignorable (vz) direction. The electron beam speed co- 
incides with the DAW seed. We gather from Fig. 3(a) 
that accelerated beam electrons have speeds 0.25c < 
\vx,e\ < 0.35c and from Figs. 4(a), 5(a) and 6(a) we have 
0.25c < \v-j:^e\ < 0.45c. We note from Table U that phase 
speed of DAWs is 0.2c— 0.26c. Thus electron acceleration 
seems to be due to the wave-particle interaction effects. 
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FIG. 8: As in Fig. 3 but for the numerical run ERI61. 
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FIG. 9: As in Fig. 3 but for the numerical run L73. 



Namely, due to the Landau damping in which wave gives 
up energy whilst accelerating particles. Particles which 
are faster than the wave give up their energy to the wave, 
while particles which are slower gain energy from the 
wave. For a given phase speed of the wave, VA,ph, in the 
Maxwellian distribution there are always more particles 
with V < VA,ph than particles with v > VA,ph- Thus, 
the net result of having Maxwellian distribution is that 
the wave gives up its momentum to the particles with 
V < VA,ph and accelerates them. As a result the accel- 



erated particle number increases around v sa VA,ph, i-e. 
a beam forms with v ~ VA,ph- It is interesting to note 
that similar conclusion has been drawn from theoretical 
calculations (29|]. It should be stressed that their calcu- 
lations were performed in the kinetic regime. Whereas 
the above conclusions are in the inertial regime (because 
of the artificially lowered mass ratio). We also note that 
similar conclusions still hold in the kinetic regime in our 
simulations (see below). Thus electron acceleration can 
be understood by simple Landau damping of DAWs de- 
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spite of the regime (inertial or kinetic). Importance of 
Landau damping in this context was also pointed out by 
Ref.fii. 

Figs. 7 and 8 correspond to runs ELI61 and ERI61. We 
see from these figures that the behaviour is similar to run 
L16. L-and R-elliptically polarised DAW with electric 
field in the non-ignorable (y) transverse direction being 
l/6th that of in the ignorable (z) direction produce weak 
parallel electron beams and weak transverse ion heating 
(no pronounced ion beams). 

Fig. 9 shows the results for the run L73 which brings 
us into the kinetic regime. By comparing this run to L16 
which is in the inertial regime we note following differ- 
ences: (i) bumps in the parallel electron velocity distribu- 
tion function, which represent the electron beam created 
by the parallel electric field, now shift to the lower veloc- 
ities which correspond to the phase speed in the DAW. 
Table H] provides DAW phase speed of « 0.1c, which in- 
deed corresponds to the bump in the f{vx,e) (in Fig. 9(a)) 
at i = 75a;~^. (ii) AH\]^f,{t) index now attains 29% which 
is about 50% larger than that of the run L16. This can be 
explained by the fact that at lower velocities there more 
particles with lower velocity. At present rrii/mf, = 1836 
is not accessible due to high computational costs. But 
from run L73 it is clear that more massive ions would 
yield even higher fraction of accelerated electrons with 
speeds comparable to the DAW speed. 

Fig. 10 is analogous to Fig. 1 but here mass ratio is 
73.44, as it corresponds to the run R73. Note that as 
we have increased ion mass from 16 to 73.44 we had to 
about double the domain length in order to prevent the 
two DAWs from colliding with each other, so the results 
are not affected by the wave interference. As can be seen 
from Fig. 10 phase mixing of DAWs develops in the ki- 
netic regime in a similar way as in the the inertial regime. 
We note that in the kinetic regime parallel currents are 
stronger 1.3 (Fig. 10(g)) compared to the inertial regime 
0.83 (Fig. 1(g)). 

We also gather from Figs. 11(a) and 11(b) that in 
the kinetic regime kink-like oscillation amplitudes in the 
electron and ion density are larger than in the inertial 
regime (cf. Figs. 1(a) and 1(b)). Also, electron ion 
heating affects larger areas around the density gradients 
(Figs. 11(c) and 11(d)). 

Figs. 12, 13, 14 (kinetic regime) are analogous to 
Figs. 4, 5, 6 (inertial regime) and correspond to the runs 
R73, EL73 and ER73. This means we can directly com- 
pare differences in the inertial and kinetic regimes. The 
conclusions that can be drawn are that in the kinetic 
regime, R-circularly polarised, L-and R-elliptically po- 
larised DAW with electric field in the non-ignorable [y) 
transverse direction exceeding 6 times that of in the ig- 
norable (z) direction produce more prominent parallel 
electron beams (^iJj|_e(t) « 35%) and transverse ion 
beams in the ignorable (vz) direction, compared both 
to the corresponding runs in the inertial regime and L- 
circularly polarised wave in the kinetic regime. Again, 
the electron beam speed coincides with the DAW speed. 



Therefore the electron acceleration can again be under- 
stood by a Landau damping of DAWs. 

Figs. 15 and 16 correspond to runs EL73i and ER73i. 
We see from these figures that the behaviour is similar 
to run L73. L-and R-elliptically polarised DAW with 
electric field in the non-ignorable {y) transverse direction 
being l/6th that of in the ignorable (z) direction produce 
weak parallel electron beams and weak transverse ion 
heating (no pronounced ion beams). 

We finish presentation of the results by considering one 
long time run L16Long which is similar to L16 but for 4 
times larger t^nd — 300a;^ . Accordingly we had to in- 
crease the domain length in x-direction 4 times. Thus in 
L16Long run we have 20000 x 200 grid size with 8 x 10^ 
particles in the simulation. The purpose of run L16Long 
is to check whether numerical simulation results pre- 
sented in runs L16, R16, EL16, ER16, ELI61 and ERI61 
have attained their respective "saturated" values, as for 
example from Fig. 3(g) or Fig. 3(h) it is not entirely clear 
what happens to the amplitude of the generated parallel 
electric field and fraction of accelerated particles in the 
density inhomogeneity regions for asymptotically large 
times. One has to realise that tend = Tbuj^i is only 
2.127 X 10~^ s (for mass ratio mi/me)- Which is few 
orders of magnitudes smaller that impulsive stage of the 
flare which is believed to be generating the considered 
DAWs. 

In Fig. 17 we plot time snapshots of various physical 
quantities at tend — 300a;^^. Fig. 17 is analogous to 
Fig. 1 but for longer times run. Commensurately we see 
larger number of full wavelength (namely 13) develop- 
ing in the driven DAW front. Again the simulation is 
stopped before the DAWs collide, which Fig. 17 clearly 
demonstrates. In Fig. 18 we show time evolution of the 
electron and ion distribution functions (Figs.l8(a)-(f)), 
amplitude of the generated E\\ (Fig. 18(g)), and AH\\^e{i) 
and AH^__i{t) indexes (Fig. 18(h)). By comparison with 
Fig. 3, we gather from Fig. 18 that parallel electron beam 
becomes more intense (Figs. 18(a)), and minor transverse 
electron heating occurs too (Figs. 18(b),(c)). Ion beams 
in the transverse direction form too but these are very 
weak. It is interesting to note that electric field sat- 
uration level did not change much (still ~ 0.03) both 
yli/ji e(t) and AH±^i{t) indexes attain the same value of 
30%. Remeber that in case of L16 AH\\e{t) attained only 
20% by the time tend = 750;" "'^. Thus we conclude that 
in run L16 the physical quantities have nearly reached 
their asymptotic values. 



IV. CONCLUSIONS 

X-ray observations of solar fiares indicate that 10"^^ — 
10*^^ per second accelerated electrons are required to pro- 
duce the observed hard X-ray emission. If the particle ac- 
celeration takes place at the loop apex, this implies that 
for the acceleration volume of 1 — 10 Mm'^ with the num- 
ber density of n = 10^® m^^^ all 100% of electrons need 



13 




D 2on *ao &no aoa iQoa 
^ : I 

-2.0=-n2 -a.0=-03 I.7C-D3 1.5c- 

E. (d) 




20D 400 6D0 SW 1000 



H 



-3,4*-n2 -1.-(a-02 7.0P-Q3 2J9-0! 
20 



%A'^ 



/ > 



o E'- ■■■ 



3DD 400 anO SKI iQoa 











B 


(b) 






aj 










1 


'■ 




15 










1 






10 










1 






5 










1 

















i 









a.W 


15 


1/ ^^\1 


ID 
S 


lUJfWiJ 


a 


L M. M .M ■■■ ■■■ .B 



20D *ao eoo aon inon 



20D 4D0 BOO BOO 1 000 
» [=/"p.l 




200 400 BOO aOQ 1D0D 

^ ■ I 

-3,Jp-0E -1.£a-0! 9.05-03 J.Da-Oa 




200 4D0 BOO BOD IQOD 

^ I 

-1,0p-0£ -J.2a-03 3.aa-03 1.1a-0a 




20D 400 eOC SOD 1 DOD 



200 4D0 BOO BOD 1 DOD 
» [=/"p.l 



-I.ie4-D0 -d.SE-OI 5.SE-D1 



-a.BE-02 I.Ie-01 



-5.4E-D1 -1.2E-01 I.Oe-01 



FIG. 10: As in Fig. 1 but for the numerical run R73. 
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FIG. 11: As in Fig. 2 but for the numerical run R73. 



to be accelerated. Given the implausibility of the latter 
acceleration efficiency, if the flare energy release triggers 
DAWs, these can rush down the solar coronal loop to- 
wards the footpoints and accelerate electrons along the 
way to chromosphere |35| . In this context and in the light 
of our initial study Ref . [221 ' ^^^^ ^^ study particle accel- 
eration by the low frequency (w = 0.3ujci) DAWs, which 
in fact are very similar to usual shear Alfven waves, as far 
as their dispersion properties are concerned. The DAWs 



propagate in the transversely inhomogeneous plasma. In 
which density is increased in the middle, across the mag- 
netic field, by factor of four and temperature is lowered 
accordingly to keep the structure in pressure balance. 
Such system mimics a solar coronal loop. DAWs generate 
parallel electric field only in the density inhomogeneity 
regions (loop edges). This parallel electric field behaviour 
coincides with the parallel current J — e{niVi — rieVe). 
The latter indicates that charge separation (due to differ- 
ent electron and ion inertia) is the cause of the parallel 
electric field generation. £'|| then accelerates electrons 
(creates electron beams) along the magnetic field, while 
E± mostly heats ions in the transverse direction. In this 
study we focus on the effect of the wave polarisation, L- 
and R- circular and elliptical, in the different regimes in- 
ertial (/3 < m^/mi) and kinetic (/3 > me /mi). Plasma 
beta is fixed at 0.02. Hence we vary the regimes by set- 
ting either m,f./mi = 1/16 = 0.0625 > /3 (inertial regime) 
and me/m,i — 1/73.44 = 0.0136 < /3 (kinetic regime). 
The mass ratio of 73.44 is not ideal (it corresponds to 
l/25th or the real mass ratio of 1836) but it is on the 
limit of available computational resources as a typical 
with m,i/m.e = 73.44 run takes 48 hours on 256 proces- 
sors (Here we present the results from total of 13 runs). 
Our main goal was to study the latter effects on the ef- 
ficiency of particle acceleration and the parallel electric 
field generation. As a result we have established: (i) 
In the inertial regime, fraction of accelerated electrons 
(along the magnetic field), in the density gradient re- 
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FIG. 12: As in Fig. 3 but for the numerical run R73. 
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FIG. 13: As in Fig. 3 but for the numerical run EL73. 



gions is ~ 20% by the time of IC driving of tend — ^So;^^ 
which allows for the 3 full wavelength to develop and is 
increasing to ~ 30% by the time of t^nd = 300aj~^ that 
allows for the 13 wavelength to develop. In all considered 
cases (see Table |ll|) ions are heated in the transverse to 
the magnetic field direction and fraction of heated ions 
(with velocities above initial average thermal speed) is 
« 35%. (ii) While keeping the power of injected DAWs 
the same in all considered numerical simulation runs, in 
the case of right circular, left and right elliptical polari- 



sation DAWs with the electric field in the non-ignorable 
{y) transverse direction exceeding several times that of 
in the ignorable (z) direction produce more pronounced 
parallel electron beams (with larger maximal electron ve- 
locities) and transverse ion beams in the ignorable (w^) 
direction. In the inertial regime such polarisations yield 
the fraction of accelerated electrons about 20%. In the 
kinetic regime this increases to 35%. (iii) In the case of 
left circular, left and right elliptical polarisation DAWs 
with the electric field in the non-ignorable (y) transverse 
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FIG. 14: As in Fig. 3 but for the numerical run ER73. 
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FIG. 15: As in Fig. 3 but for the numerical run EL73i. 



direction several times less than that of in the ignor- 
able (z) direction produce less energetic parallel electron 
beams. In the inertial regime such polarisations yield 
the fraction of accelerated electrons about 20%. In the 
kinetic regime this increases to 30%. (iv) The parallel 
electric field that is generated in the density inhomo- 
geneity regions is independent of the regime (whether 
kinetic or inertial), i.e. independent of electron- ion mass 
ratio and stays of the order 0.03ujpeCme/e which for con- 
sidered solar flaring plasma parameters (n = 10^^ m~^ 



and T = 6 X lO^K at the domain edges) exceeds Dreicer 
electric field Ed = nge^ lnA/(87re§fcBr) (with Coulomb 
logarithm InA = 17.75) by a factor of 6.4 x 10*. Al- 
though our PIC simulations do not include particle col- 
lisions, and thus Dreicer electric field is not directly rel- 
evant because this concept is essentially collisional, the 
fact that the generated parallel electric field exceeds Dre- 
icer value by eight orders of magnitude indicates that the 
accelerated electrons would be in the run-away regime. 
The fact that in the present simulations Eu amplitude is 
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FIG. 16: As in Fig. 3 but for the numerical run ER73i. 
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FIG. 17: As in Fig. 1 but for the numerical run L16Long. 



independent of nii/me contradicts with the scahng law 
(i?ll ex l/{mi/me)) found by Ref. 3(^]. The difference can 
possibly explained by the fact that Ref. 36] considered 
cold plasma approximation i.e. all particles had the same 
velocity (with the distribution function being a delta- 



function rather than a Maxwellian). The latter would 
have prevented the possibility of having wave-particle in- 
teractions, which as found in this study are rather im- 
portant, (v) Electron beam velocity has the phase ve- 
locity of the DAW. This can be understood by Landau 
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FIG. 18: As in Fig. 3 but for the numerical run L16Long. 



damping of DAWs. DAWs give up their energy to elec- 
trons via collisionless Landau damping. Because for the 
DAWs with ui = O.Swci their phase seed is approximately 
the Alfven speed, this, combined with the suggestion by 
Ref.[35l| that Alfven speed can be as high as few xO.lc, 
for the realistic mass ratio of rrii/me = 1836, would 
readily provide electrons with few tens of keV. e.g. for 
Va = 0.3c the electron energy would be 23 keV. (vi) As 
we increased the mass ratio from mi/me = 16 to 73.44 
the fraction of accelerated electrons has increased from 
20% to 30-35% (depending on DAW polarisation). This 
is because the velocity of the beam has shifted to lower 
velocity (compare e.g. Fig. 3(a) to 9(a)). Since there are 
always more electrons with a smaller velocity than higher 
velocity in the Maxwellian distribution, for the mass ratio 
mi/me = 1836 the fraction of accelerated electrons would 
be even higher than 35%. (vii) DAWs generate significant 
density and temperature perturbations that are located 
in the density gradient regions (solar coronal loop edges). 
These perturbations of density and temperature should 
in principle be detectable in future observations because 
current time cadence of space instruments is well above 



The overall conclusion is that DAWs, propagating in 
the transversely inhomogeneous plasma, in the Alfvenic 
(w < LUci) range of the plasma dispersion can effectively 
accelerate electrons along the magnetic field and heat the 
ions in the transverse direction. This is bound to have 
applications in many fields - Earth magnetosphere and 
laboratory plasma, in addition to the solar flare context 
emphasised here. 
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